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Abstract 

Motivated by recent research of aging in E. coli, we explore the effects of aging 
on bacterial fitness. The disposable soma theory of aging was developed to explain 
how differences in lifespans and aging rates could be linked to life history trade-offs. 
Although generally applied for multicellular organisms, it is also useful for exploring life 
history strategies of single celled organisms such as bacteria. Starting from the Euler- 
Lotka equation, we propose a mathematical model to explore how a finite lifespan 
effects fitness of bacteria. We find that that there is surprisingly little loss of fitness 
when the bacterium has limited opportunities to reproduce. Instead, the fitness gained 
each time the bacteria reproduces decreases rapidly, and investing resources to survive 
to reproduce the first few times is likely more advantageous than investing additional 
resources to try to maintain cell integrity longer. 
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1 Introduction 



Living organisms exhibit an astonishing variety of hfespans and hfe histories. A giant sequoia 
tree can hve for thousands of years, producing only a few offspring over its hfetime, while 
a dandilion survives only for a season, but produces thousands of seeds, often resulting in 
hundreds of new plants the next year. Among a single genus, such as rockfish (Sebastes)^ 
species exhibit lifespans ranging anywhere from 10 to 200 years dLove et al.,l002| ). Even 
bacteria appear to age. New research by Stewart et al. ( [Stewart et al., 2005D indicates that 
even Escherichia coli which appear to divide symmetrically actually divide into one "old" 
and one "young" cell. Current theories of aging seek to combine principles of evolution with 
theories from microbiology about how damage that accumulates in a cell makes it unable to 
function. One such theory is called the disposable soma theory of aging ( |Kirkwood,T98Tl 
[Finch aiiH" Kirkwood, 2000l |D renos and Kirkwood, 200*51 ). This model predicts that an or- 
ganism will allocate resources for reproduction and maintainace of the soma in such a way 
that maximizes the dispersal of its genes. We expect that the optimal allocation strategy 
will not be one that allows an organism to maintain the soma indefinitely. This is because 
organisms have a finite amount of energy to use for all life functions; if they use all of the 
energy for maintenance of the soma, then there might not be enough energy to reproduce, 
and vice versa. This theory also predicts how natural mortality is expected to infiuence the 
life span and reproductive schedule of an organism. For instance, if an organism experiences 
high natural mortality, then its resources are better invested in offspring than soma. Higher 
natural mortality also results earlier maturation, so they will be less likely to die before hav- 
ing an opportunity to reproduce. Organisms with lower natural mortality instead maintain 
the soma for longer, produce offspring less frequently, and experience longer lives. 
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2 Fitness Models of Bacteria 



Single celled organisms, or other organisms where the soma and the germ line are not sep- 
arate, were once expected to be immortal ( Williams, 1957D . The distinction between or- 



ganisms that age and those that do not has more recently been hypothesized to depend 
upon asymmetry in reproduction ( [Partridge and Barton, 1993| ). Recent research by Stewart 
et al. dStewart et al., 2005D indicates that even bacteria that appear to divide symmetri- 
cally, such as E. Coli^ actually produce functionally asymmetric cells during cell division. 
They consider one of the cells to be an aging parent cell that produces offspring that are 
''rejuvenated" and find evidence evidence that these older cells reproduce more and more 
slowly as they age, and may even stop reproducing after a certain number of divisions 
dStewart et al., 2005D . Quantitative models can be useful for exploring how evolutionary 
trade-ofl[s shape aging and senescence in these simple organisms. Here we explore the effect 
of finite lifespan on bacterial fitness when doubling rate and environment are constant. 

We explore the effects of life history choices on our selected fitness measure, the intrinsic 
rate of natural increase, denoted by r, in populations with age dependent reproduction and 
mortality schedules, living in a constant environment using the Euler-Lotka equation: 



oo 



e ''''Ixbxdx. (1) 







Here the probability of surviving to age x is denoted as Ix and the number of offspring born 
to an individual of age x is bx- 

Kirkwood ( [Kirkwood, 1981D first proposed a simple model of bacterial fitness for cells 
that divide perfectly symmetrically, based upon (P) , with appropriate choices of bx and Ix for 
a clonally reproducing population. One simple choice for Ix is an exponentially decreasing 
survival probability, Ix = e~^^, where m is the mortality rate. The birth rate will be related 
to the doubling time, T. If an individual bacteria survives to time T, it divides, and two 



identical offspring result. At the end of this time period the original bacteria is essentially 
"dead" and the daughter cells remain. An appropriate "birth" rate would therefore be 
bx = 2(5(x — T), where 5{x — T) is the Dirac delta function-^. With these expressions for Ix 
and bx (P) becomes: 

2e-(m+r)T _ ^ (2) 

Figure shows this functional relationship between r, m and T graphically. As m and T 
increase, r decreases. As T ^ 0, r ^ oc, regardless of the value of m. As m and T increase, 
r decreases. 




Figure 1: The intrinsic rate of natural increase, r, as a function of the doubling time, T, and the mortality, 
m. 

We denote the fraction of resources allocated for growth and reproduction by p, and the 

fraction alloted for maintenance and survival by 1— p. Following Kirkwood ( |Kirkwood, 1981D , 
^The Dirac delta function is defined as a unit impulse at some point xq such that: 

6{x — xq) =0, X xo 




S{x — xo)dx — 1, 



and given an arbitrary function f{x): 

/oo 
f{x)S{x - xo)dx = f{xo). 
-co 
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we can parameterize the mortality, m, and doubling time, T, in terms of p as: 

r(p) = ^ (3) 

m{p) = (4) 

Here, Tq can be thought of as the minimum possible time it would take for the bacteria to 
reproduce if it allocated all of it's resources to growth; mo is the minimum mortality of the 
bacteria if all resources are allocated for survival. Solving for r in (|2j) with the expressions 
for T and m in (jHlEj) yields: 

r = — ln2- -. 5 

To 1 - p 

Maximizing © with respect to p gives the optimal resource allocation, p*: 

and the corresponding value of r^axi 

Tmax = ^ (ln2 - 2(moToln2)^) . (7) 

If To and mo are constants, r{p) will follow curves similar to those depicted in Figure [21 
We can see from Eqn. (jSj) and Figure (21 that the optimal strategy will never be to dedicate 
all resources to reproduction, i.e. p* < 1. In an environment with low mortality, it can 
be optimal to invest most resources in reproduction, as shown in Figure Et^, even if the 
generation time is relatively long (lowest curve). However, if the mortality is high and/or 
the doubling time is long, the optimal resource allocation may be to use more resources for 
survival (Figure Eb). 



6 




1" p 



(a) 



r(p) 

5 \ 



-2.5 
-5 
-7.5 
-10 




To=0 


06 


To=0 


1 


To=0 


25 


To=2 





(b) 



Figure 2: r as a function of p for four values of Tq = (2,0.25,0.1,0.06) for (a) mo = 0.015 (b) mo = 0.9. 
High values of To correspond to the lowest curves in each plot, and small To to higher curves. 

3 The Implications of Bacterial Aging 

We now explore how the intrinsic rate of natural increase for the aging bacteria, denoted by 
f, is impacted by a finite lifespan. Starting from the Euler-Lotka Eqn. (P) we first modify 
the birth rate, h^^ to take into account a functional asymmetry in cell division. The cell is 
able to divide and produce a single oflFspring in a given fixed doubling time T, then can live 
to divide again later. We define cellular "age", a = 1, 2, ... , amax-> as the number of times 
the cell has doubled, where a^^x is the maximum number of times it can split (in analogy 
with the Hayfiick limit ( |Arking, 19981 |Finch and Kirkwood, 2000| )). The birth rate is then 
a sum of delta functions spaced a distance T apart: 



hx 5{x — aT). 



(8) 



a=l 
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This, together with the previous expression for inserted into (P) gives: 

poo dmax 

1= / e-^^e-*"^ V - ar)dx 

a=l 

a=l 

Evaluating the sum in Eqn. dHj) results in 

^{r+m)T — 2 — ^-(^^^)Tamax ^ j^^Q) 

This equation has two solutions for a^ax ^ 1- The trivial solution exists when f = —m. 
All other solutions depend upon the maximum age, a^ax^ the mortality rate, m, and the 
doubling time, T (see Figure Ej). The value of f varies considerably depending upon the 
combination of these three parameters. Variation of T and a^nax seem to have the most 
impact upon f, as shown in Figure Ett, whereas variations in m act to shift f up or down 
a fixed amount when T is held constant (Figure Eb)- Variation in a^ax has a much larger 
impact upon r when the value of a^ax is small. When a^ax co, Eqn. (^Uj) reduces to Eqn. 
(|2I), and f ^ r. In other words, infinite reproductive lifespan and perfectly symmetrical cell 
division are equivalent. 

Since as ajnax increases f ^ r, we can approximate the second solution of Eqn. ffTUj) in 
terms of r for large values of a^ax using Newton's method: 

2 ^rnax 

r ^ r — —- . (11) 

Tln2 + a^ax2--max ^ ^ 



As can be seen in Figure this approximation is remarkably good, over a wide range of 
values of ajnax • It also gives us some idea of the values of a^nax that are "large" . Where the 
approximation is valid, a^ax must be large enough to make the difference in fitness between 
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Figure 3: f vs amax for (a) m = 0.5, T = (0.25, 0.5, 0.75, 1) and (b) T = 0.5 and m = (0.5, 0.25, 0.1, 0.01) (c) 
Effects of amax = (2, 3, 5, 10) on f(T) for m = 0.9 

the aging bacteria and the immortal bacteria, r — very smaU. 

Since the fitness of the aging and immortal bacteria are very close when a^ax is large, 
we would expect that the allocation strategies would be similar as well. Let us denote the 
allocation strategy of an aging bacteria by p. Substituting dHj) and into ffTUj) gives a 
relationship between f, p, mo, Tq, and ajnax' 

exp ( (f + ) = 2 - exp ( -(f + -^^)5a^«^ ). (12) 

V ^- P P J \ ^- P P J 

Unlike the baseline model without aging, Eqn. ffT2j) does not have an exact closed form 
analytic solution. Figure El shows a numerical solution for the fitness of the aging bacteria, 
f, as a function of fraction of resources allocated for reproduction, p, for various values of 
<^max5^0 5 and Tq. The values of f increase drastically between the lowest values of a^ax- 
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Figure 4: (a) Analytic approximation (dashed lines) and numerical solution (solid lines) of r vs amax foi" 
m = 0.5 and T = (0.1,0.25,0.5,0.75) (b) The difference between the numerical solution and analytic ap- 



proximation, f 



for the same values as in (a) 



However, for higher values of a^^x? ^ hardly varies and there is a lower return in the fitness 
for more opportunities to reproduce. Additionally, notice that as a^ax increases, the peak in 
f becomes more pronounced, and shifts towards higher values of p. Thus, bacteria that can 
double many times gain more fitness from focusing resources for reproduction than those 
that cannot. 

We can quantify how aging shifts the optimal allocation strategy away from the value 
determined in Eqn. © for immortal bacteria. If we differentiate Eqn. ([T^ with respect to p, 
and evaluate it at the optimal strategy, p*, so that ^^^^ ^ = 0, we find an implicit expression 
for the maximum fitness, f(p*) = Trnax^ as a function of p*: 

((2p* - l)mo - fmax{l - p"f) [e(^— +^(^*))^(^*) + a^,,e-(^— +"^(^*»^(^*)"— 



]=0 (13) 
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Figure 5: Effects of amax = (1,2,5,10,20) on r{p) and r{p) (solid black line) for (a) mo = 0.7, Tq = 0.45; 
(b) mo = 0.7, Tq = 0.25; (c) mo = 0.9, To = 0.1; (d) mo = 0.9, To = 0.01; (e) mo = 0.01, To = 1; and (f) 
mo = 0.5, To = 0.1. 

Since the sum of exponential terms in f[T!?j) can never be equal to zero, we conclude that 



(2p*-l)mo- w(l-p*)^ = 



and therefore r^riax is related to p* by 



. i2p* - 1) 



(14) 



It is straightforward to confirm that this expression holds for the case when amax oo. 
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Figure 6: Close-ups of the maxima of r{p) for amax = (5, 10, 20) and r{p) (solid black line) for (a) mo = 0.7, 
To = 0.25 (corresponding to Figure [SK ) (b) mo = 0.9, To = 0.01; (corresponding to Figure [Sli ) 

As previously discussed, as the maximum cellular age, a^ax^ increases, the difference between 
the fitness of the immortal bacteria, r, and the fitness of the aging bacteria, f, decreases. In 
particular, for large values of a^ax we saw that r ^ r. We can also see in Figure El that when 
(^max ^ 10 the peaks in the fitness curves are very close, so that r^nax ^ ^max^ and occur at 
about the same value of the allocation strategy, i.e. p" . We would like to quantify the 
magnitude of the change in optimal allocation strategy as a function of amax- We begin by 
looking at the difference between Vmax and fmax'- 

Tmax - fmax = r{p'') - r{p') + r(p*) - f(p*). (15) 
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Writing fjl5j] . in this way allows us to use our approximation from Eqn. fjllj) to determine 
the value of a small parameter, e, in our approximation. Substituting Eqns. f[TT|) . and 
ffTlj) into ([TK|) and rearranging gives: 



(2p* - 1) 

(l-p*)2 



-mo f» 



p* In 2 
To 




mo , p* 



where e = 



In 2+a; 




. Solving this expression for p*, we find: 



1 - 



moro y 

ln2-ey 



1 / ln2 

2 V"^oro 



(16) 



where 0(6^) denotes terms that are of order or smaller. From this we also have an 
approximate expression for Trnax^ similar to ([Zj): 



This analytic approximation gives results that are very close to those obtained numer- 
ically. For instance, for Tq = 0.05 and mg = 0.5, the difference between the analytic 
approximation and numerical solution when a^^x = 5 is fairly large, as we might have 
guessed since this is when r — f is largest (Figure Hj). In this case the numerical solu- 
tion is {p^^fjnax) = (0.8077,8.31959) compared to the analytic approximation {p^^fmax) = 
(0.804836,8.00324). As amax increases, the approximation improves. For instance, when 
dmax = 20, the numerical solution is (p*, Trnax) = (0.8101, 8.58945), and the analytic approx- 
imation is (p^^rmax) = (0.810086,8.59738), a difference of less then 0.1% in r^^ax and less 
than 0.002% in p*. 




(17) 
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4 Discussion 



Studying the impacts of aging and resource allocation on fitness for bacteria is appealing 
for a number of reasons. Bacterial systems are relatively simple. They also are ideal for 
experimental manipulation. Metabolism and repair activity should be fairly straightforward 
to measure, giving indications of how bacteria allocate resources. Genetic manipulation of 
bacterial systems is also possible, which allows more direct measurement of the parameters 
in the model. Because generation times are short, bacterial systems also allow for replicate 
experiments at a reasonable cost. Thus, bacterial systems are ideal model systems for ecology 
JJessup et al., 2004D . 

Here we have proposed a simple mathematical model to explore how aging, in the form of 
finite life span, effects fitness in bacteria. This allows us to explore, explicitly, predictions of 
the disposable soma theory of aging for a system that is fairly simple and easy to manipulate. 

This model provides valuable insight into which trade-oflFs have the most impact on 
bacterial fitness. Given the assumptions about the form of Ix in the Euler-Lotka equation, as 
well as the relationships between mortality, doubling rate, and resource allocation, we found 
that the ability to manipulate the doubling time T has the greatest impact upon bacterial 
fitness. If the doubling time is short, even if the mortality is high, as is shown in Figure 
Ell, the bacteria's fitness is considerably higher than in a system with lower mortality, but 
longer doubling time (Figure andEt). Thus we expect for there to be strong selection for 
lower values of the minimum doubling time, Tq. 

We also find that that a bacteria experiences surprisingly little loss of fitness when it 
has a finite number of opportunities to reproduce (Figures El and E|). For most combinations 
of the minimum mortality rate, mo, and minimum doubling time, Tq, a^nax ^ 5 or 10 is 
large enough to confer almost exactly the same amount of fitness to the bacteria as an 
ability to reproduce indefinitely (Figure E|)- Most of the bacteria's fitness is gained the first 
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few times the it doubles, and the amount gained in each subsequent doubhng decreases 
rapidly. Therefore we expect that if there were a resource cost associated with increasing 
the maximum number of possible doublings then investing resources to survive to double the 
first few times would be better than investing additional resources to try to maintain cell 
integrity indefinitely. 

This simple model also has the advantage of being analytically tractable. Although 
numerical results for complex models are fairly easy to obtain it can be difficult to understand 
and quantify the roles of model parameters. Good analytic approximations, such as those 
found above for r^Tmax^ P*, allow us to explore the effects of parameter variation in a 
very concrete way that might not be available for more complex models. 

This model is only a first step in understanding how aging impacts bacterial fitness. 
Other mechanisms are likely important in determining the optimal life history strategy. The 
model does not consider the effects of variable environmental conditions, density effects, or 
of increasing doubling times as a function of age. Also, here we assume no dependence of 
amax on p, so the trade-offs between mortality, reproduction, and lifespan are not explicitly 
explored. These could be important factors in the development of bacterial life histories. 
However, in spite of these limitations, this model gives insight into why bacteria do not 
exhibit infinite lifespans, and suggests directions for further exploration. 
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